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ABSTRACT 

Context. Water vapor emission at 22 GHz from masers associated with star-forming regions is highly variable. 
Aims. We present a database of up to 20 years of monitoring of a sample of 43 masers within star-forming regions. The 
sample covers a large range of luminosities of the associated IRAS source and is representative of the entire population 
of H2O masers of this type. The database forms a good starting point for any further study of H2O maser variability. 
Methods. The observations were obtained with the Medicina 32-m radiotelescope, at a rate of 4-5 observations per 
year. 

Results. To provide a database that can be easily accessed through the web, we give for each source: plots of the 
calibrated spectra, the velocity-time-flux density plot, the light curve of the integrated flux, the lower and upper 
envelopes of the maser emission, the mean spectrum, and the rate of the maser occurrence as a function of velocity. 
Figures for just one source are given in the text for representative purposes. Figures for all the sources are given in 
electronic form in the on-line appendix. A discussion of the main properties of the H2O variability in our sample will 
be presented in a forthcoming paper. 
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1. Introduction 

Since the discovery of 22 GHz H2O maser emission associ- 
ated with young stellar objects (YSOs) within star-forming 

■ regions (SFRs), variability of maser emission is well-known. 
. Changes as large as several orders of magnitude in the 

■ maser emission have been observed (e.g., Little et al. 119771 

■ Liljestrom et al. 119891 Claussen et al. 11996! Comoretto et 
. al., 1990, Wouterloot et al. ll995p . At the same time, velocity 
' drifts of individual components of up to a few km per 

year have also been reported (e.g.. Brand et al. l2003p . The 
variability can be slow or burst-like and covers all ranges of 
timescales, from hours-days to months-years. In the present 
study we are mainly concerned with the latter and, for the 
first time, deal with a large sample of sources (43) and with 
a time-span of up to 20 years, with 4-5 spectra per year. 

Due to the large amount of telescope time required 
to follow the evolution of H2O maser emission, inevitably 
the variability is more easily monitored through single 
dish observations. In fact, besides our campaign only the 
Pushchino 22-m single-dish maser patrol covers a compa- 
rably long period (e.g., Rudnitskij et al. 2007). Furuya et 
al. (|2003p carried out a multiepoch 22 GHz H2O maser 
survey towards 173 low-mass YSOs (Class to Class HI 
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sources) using the Nobeyama 45 m telescope. This was the 
first complete water maser survey towards Class sources in 
the northern sky. However, their observations extend non- 
uniformly over a period of only three years. 

It is well-known that the variability depends strongly on 
the luminosity of the SFR (usually derived from the associ- 
ated IRAS source). For instance, there is a minimum lumi- 
nosity (~25 Lq) below which H2O maser emission may be 
present only for about one third of the entire duration of the 
maser activity (Persi et al. il994i Claussen et al. J.996). The 
results obtained so far suggest that the lower the SFR lumi- 
nosity, the higher the observed degree of variability of the 
H2O maser emission. Higher luminosity sources may show 
more steady components. Consequently, a large sample of 
sources is needed to better understand the dependence of 
the variability on other parameters of the SFRs, in par- 
ticular on their luminosities which cover a range of several 
orders of magnitude, and on their (molecular) environment. 

The use of a single dish instrument is of course a limita- 
tion because interferometric observations clearly show that 
H2O masers within a SFR often consist of many spatially 
separated, unresolved components, generally clustered in 
groups (usually with different velocities), which in most 
cases cannot be separated in single dish observations (e.g., 
Forster & Caswell[TnHni [TO^ Tofani et aLlMF). However, 
Felli et al. (2006J reported a case in which single dish ob- 
servations were able to separate the evolution of spatially 
distinct components. Nevertheless, frequent single dish ob- 
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servations have the advantage of being able to foUow poten- 
tial velocity drifts of individual velocity components, which 
can be easily identified in the spectra, and thus give a bet- 
ter view of the dynamics (e.g., accelerations) occurring in 
the circumstellar environment of the YSO. 

In two earlier papers (Valdettaro et al. 120021 Brand et 
al. [2003) , we set the basis for a systematic study of the H2O 
variability in SFRs with the Medicina 32-m radiotelescope 
using a sample of 14 sources which covered a wide range of 
luminosities and had been observed for more than 10 years. 
Here we report on a larger sample of 43 SFRs (including 
the former 14 sources) and increase the time coverage up 
to 20 years. The database is presented in the form of plots 
of the calibrated spectra and additional plots of derived 
quantities, in a form that can be easily accessed through the 
web. In a forthcoming paper, we will analyze the properties 
of II2O variability of our sample. 



2. The sample 

A large sample of sources selected from the Arcetri Atlas 
(C esaron i et al. HMEi Comoretto et al. fTM)! Palagi et 
al. [T9931 Brand et al. [TOMl Valdettaro et al. [2001]) has 
been monitored regularly since the beginning of the op- 
erating period of the Arcetri digital autocorrelator at the 
Medicina 32-m radiotelescope (1987), at a rate of 4-5 ob- 
servations per year. Our goal was to study the dependence 
of the variability of H2O masers in SFRs on the luminosity 
of the associated YSO, assumed to remain constant during 
our patrol. It should be pointed out that the FIR luminos- 
ity is derived from IRAS observations and, consequently, is 
representative of the entire luminosity of the SFR and not 
necessarily of the YSO that powers the II2O maser. In fact a 
SFR may have a very complex structure, with several YSOs 
in various evolutionary stages and with different luminosi- 
ties present within the same SFR, as for instance in W75-N 
(Hunter et al. ll994p . Variability might also depend on other 
characteristics of the SFR such as the presence of molecular 
outflows (Felli et al. 1992), its IRAS colours and the pres- 
ence or absence of an UC Hii region (Palla et al. il99l[|1993l 
Code lla et al. [19941 [ISMl Codella & FeUi [19951 Codella & 
Palla [ISMI). 

The present sample of 43 H2O masers was chosen out of 
the several hundred known H2O masers of SFR type (e.g., 
Comoretto et al. I1990p using the following criteria: 

(1) the maser sources were to be of the SFR- type, i.e. as- 
sociated with a SFR, to be distinguished from those as- 
sociated with late-type stars. The SFR classification of 
Palagi et al. ( fT993l ) and Valdettaro et al. (jSOOP) is based 
on IRAS colours. The location of the sources of our sam- 
ple in an IRAS colour-colour plot is shown in Fig. [1] and 
compared with that of a much larger sample of SFRs 
with H2O maser emission (a merger of Comoretto et 
al. [T9901 Brand et al. I T994I and Valdettaro et al. [200Tj) . 
The dashed box defines the region occupied by sources 
with colours characteristic of ultracompact Hii regions 
(Wood & Churchwell [19891). To double-check the true 
nature of the sources exciting the H2O masers, other 
morphological indicators of the SFRs were also searched 
for in the literature, e.g. association with dense molecu- 
lar cores (e.g Cesaroni et al. I1999p . presence of UC Hii 
regions, presence of molecular outflows and other types 



of masers. All investigations further support the SFR 
identificatiorQ. 

(2) the luminosities of the associated IRAS source was to 
cover as large a range as possible; 

(3) the sample was to be sufficiently large to be represen- 
tative of the entire population of H2O masers within 
SFRs; 

(4) the sources were to be suitably distributed in right as- 
cension, to avoid scheduling problems due to clustering 
of sources in the inner region of the Galactic plane in 
the large sample; 

(5) the peak fiux density was to be greater than a few Jy 
so that the source could be easily observed with a 5 
min on-source integration time with the Medicina 32-m 
radio telescope; 

(6) the number of selected sources was not to be too large 
to allow observations in a < 2-3 day session; 

(7) finally, in order to prevent confusion in the single beam 
from nearby, but unrelated, maser sources, we selected 
sources from the Arcetri Atlas for which no other masers 
had been reported within a circle of several arcminutes 
with single dish observation^. 

The main source parameters are given in Table [l] which 
lists: 

(1) a running number; 

(2) the source name; 

(3) the associated IRAS source, when available; 

(4) J2000.0 right ascension; 

(5) J2000.0 dechnation; 

(6) Vci, the velocity of the associated molecular cloud rela- 
tive to the LSR, derived from the literature using tracers 
sensitive to high molecular density, such as NH3 and CS 
(in km s^^) 

(7) d, the distance taken from the literature (in kpc); 

(8) LiR, the integrated FIR luminosity (in L0), usually de- 
rived from IRAS data or, otherwise, taken from the lit- 
erature; 

^ We note that in SIMBAD, AFGL2789 is identified with 
"V645 Cyg, pulsating variable star". However, as stated by 
Clarke et al. (|2006p . V645 Cyg is a "relatively unembedded 
young massive star, with a high- velocity wind and an associated 
optical and molecular outflow", and it "may represent a rela- 
tively rare class of transition objects between a genuine massive 
YSO and a normal young Oe-type star in a weak Hll region". 
Its presence in our sample is therefore justified. 

^ A possible exception might be OMC 2 and KL IRC2 which 
are separated by 13. '3. During the mega burst of KL IRC2 in 
January 1999, in which the flux density increased by at least 
a factor 10'^ and reached ~ 2.6 10^ Jy, we observed a similar 
burst in OMC 2 at the same velocity (around 8 km s~^) and 
with a flux density of ~ 3.6 10"* Jy. Beam shape mapping indi- 
cates that at this distance from the beam axis the attenuation 
should be about 10~^. Consequently, a side- lobe effect of the 
observed intensity is expected. However, apart for the period of 
the mega- burst, all the remaining data for OMC 2 should be un- 
affected by the presence of the close-by source KL IRC2, at least 
outside the velocity range close to the KL IRC2 peak emission, 
approximately for velocities less than 6 km s~^or greater than 
10 km s~^. In fact, 1) when the flux of KL IRC2 is at a level of 
2-3 10'^ Jy there are OMC 2 spectra with no signal above the 
noise level of 1-2 Jy, 2) there are features observed in OMC 2 
outside the above-mentioned velocity range which do not match 
with any emission in KL IRC 2. In conclusion, the data of OMC 
2 should be used with the caveat of possible KL IRC2 side lobe 
effects in the range of velocities with strong KL IRC2 emission. 
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(9) the date of the first observation; 
(10) references for the distance. 

The IRAS luminosity was obtained by adopting the dis- 
tance from Col. 7, assuming an emissivity proportional to 
the frequency, correcting the observed IRAS flux densi- 
ties for the derived colour temperatures between two ad- 
jacent bands, and extrapolating the fluxes to 6 and 400 fim 
(Wouterloot et al^MB- 

The sample covers a range of IR luminosities from 20 L0 
to 1.5 X 10^ L0. These values bracket the entire luminosity 
interval of the regions where H2O maser of the SFR- type 
are found (Palagi et al. 1993) . In the histogram of Fig. [2] 
the distribution of the IR luminosities of the sources in our 
sample is compared with that of a much larger sample of 
SFR with II2O maser emission (Palagi et al. [T993) . Both 
histograms are normalized to their respective total number 
of sources. The small decrease of sources in our sample at 
large luminosities is probably due to the clustering of bright 
masers in the inner region of the Galactic plane which is 
less sampled according to criterion (4). The slightly larger 
percentage of sources in our sample at low luminosities is 
due to II2O masers added after the publication of the Palagi 
et al. P^M)) sample. 

Similarly, in Fig. [3] the distribution of H2O maser in- 
tegrated fluxes of our sources in the last observing run 
(February 2007) is compared with that of a much larger 
sample of SFR (a merger of Comoretto et al. 119901 Brand 
et al. 119941 and Valdettaro et al. l200ip . Both histograms are 
normalized to their respective total number of sources. 

Figure [3] shows that in Feb. 2007 there were percentage- 
wise, more sources with high integrated flux densities in 
our present sample than there are in the comparison cata- 
logue. Likewise there is a smaller percentage of sources with 
low integrated flux densities. For a large range in Ljr the 
general trend is for the average integrated flux density to 
increase with increasing Lib, (e.g.. Brand et al. (|2003p : this 
trend will be further investigated in a forthcoming paper) 
and the difference between the distributions shown in Fig. [3] 
reflects the fact that we have biased our present sample to- 
wards sources of high Ljr, in order to have a good chance 
of yielding a maser detection in a reasonable amount of 
integration time. On the other hand, sources of lower Ljr 
usually have a lower flux density, and do not have a 100% 
detection rate (e.g.. Brand et al. I2007p . This, in combina- 
tion with our present sample's bias towards sources with 
higher Ljr can be seen as the main cause of the difference 
between the histograms in Fig. [3] for the lower integrated 
flux density bins. Furthermore, on the low flux density side, 
the difference between the histograms is enhanced by the 
way sources from the comparison catalogue were counted (if 
there was more than one observation, we always selected the 
one in which the source was detected). Note that although 
in Fig. [3] we show the distribution of integrated flux densi- 
ties for the observing session of Feb. 2007, any other session 
would have given the same statistical result. In summary 
and within the limitations inherent in the comparison, we 
confirm that our sample is fairly representative of the entire 
population of SFRs. 

Ideally, from the observed maser variability, we would 
like to infer the properties of the exciting YSO and of its 
surroundings. However, while for II2O masers from late- 
type stars the association is unambiguous (i.e. there is only 
one star pumping the maser), for SFR- type masers, sev- 
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Fig. 1. Colour-colour (Ca-f, = log(Fa/Fb), F^ = IRAS flux 
in band a fim) diagram of the IRAS counterparts of the 
sources listed in Table 1 (o), compared to that of a larger 
sample of H2O masers of SFR type (a merger of Comoretto 
et al. [T9901 Brand et al. [T9941 and Valdettaro et al. [200T|) 
( X ) . The dashed box defines the location of Ultracompact 
Hii regions (Wood & Churchwell [TOggi) 



eral distinct YSOs may simultaneously be present within 
our beam. This ambiguity can only be solved with high- 
resolution interferometric observations, but even if in many 
cases these are available, they cover at most a few epochs 
and there is no guarantee that the spatial/ velocity situa- 
tion of the maser spots remains the same over long time 
intervals, given the large variability observed. 



3. Observations 

The present study is based on observations carried out 
with the Medicina 32-m radiotelescope0 (HPBW --1.'9 at 
22 GHz). This is primarily set up to be used for VLBI mea- 
surements, and therefore the front-end is optimized for such 
work. 

The monitoring program is ongoing but the time in- 
terval considered in this work spans 20 years, from March 
1987 to February 2007, with shorter time coverage for a 
few sources which were added to the list at a later date. 
During this period, various parts of the whole system (an- 
tenna, receiver, autocorrelator) were improved, increasing 



^ The Medicina 32-m VLBI radiotelescope is operated by 
INAF-Istituto di Radioastronomia. 
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Table 1. The H2O maser sample 
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2 
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12-JUN-1987 


29 


39 
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22543+6145 


22 


56 


18 


1 


+62:01:46 
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2 


5 


10* 


ll-FEB-1987 


30 


40 


WB89-234H2O 


23004+5642 


23 


02 


31 


8 


+56:57:44 


-53.5 


5.6 


9 


6 


10^ 


17-SEP-1995 


31 


41 


Sh 2-158 


23116+6111 


23 


13 


44 


7 


+61:28:10 


-56.9 


2.5 


2 


5 


10^ 


12-JUN-1987 


32 


42 


IRAS 23139+5939 


23139+5939 


23 


16 


10 


3 


+59:55:29 


-44.0 


3.5 


1 





10* 


12-JUN-1987 


33 


43 


IRAS 23151+5912 


23151+5912 


23 


17 


20 


8 


+59:28:47 


-54.7 


3.5 
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9 


10* 


12-JUN-1987 


33 



Distances from: (1) Georgelin (fT975ll : (2) Hach isuka et al. ([2006]) : (3) Herbig & Jones (fT983)l : (4) Sneh et al. (pM]) : (5) adopted 
distance by Fukui et al. (|1986p . Meehan et al. I|1998|) quote d = 0.45 kpc; (6) adopted distance, for a recent review of distance 
determinations to the Orion region see Jeffries l|2007|l : (7) Palagi et al. I|1993|l : ( 8) Mo Unari et al. Ip002)l : (9) Evans & Blair lfT98T|) : 
(10) Anthony- Twarog (1982); (11) Shepherd fc Chu rchwcU (1996); (12) Racine (fT968)) : (13) Lada & Woode n (|T979|) ; (14) Racine & 
van den Borgh (1970J; (15) Hunter & Massey (fTMTl) : (16) M. Honma, private communication; (17) Walker (T^SS); (18) ChurchweU 
et al. ((1990); (19) kinematic distance, using the rotation curve of Brand & Blitz (1993); (20) Valdettaro ct al. (2002); (21) Honma 
et al. (PPOS); (22) Y. Xu, private communication; (23) Wilking et al. (PMS]) : (24) van dcr Tak ct al. (1999); (25) Hunter et al. 
(fT994)) : (26) C larke et al. (12006}; (27) de Zeeuw etal. (1999); (28) Crampton & Fisher (1974); (29) Molinari et al. (|T996]) ; (30) 
Blaauw et al. (|1959ll : (31) Brand & Wouterloot p998p : (32) L. Moscadelli, private communication; (33) Wouterloot & Walmsley 
(p86)l : (34) Armandroff & Herbst iflMT]) . 

the sensitivity and spectral resolution of the observations. 512-channel digital autocorrelator; in 1991 the number of 
A detailed description of the system and its improvements channels was increased to 1024. In the same year a HEMT 
can be found in Comoretto et al. (1990) and in Brand amplifier replaced the GaAs FET front-end, reducing the 
et al. (|1994p . respectively. In 1989 realigning of the an- system temperature. In 1991 an active sub-reflector control 
tenna surface resulted in an improvement of the efficiency increased the gain at low and high elevations, 
at 22 GHz from 17 to 38%. The original backend was a 
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Fig. 2. The distribution of IR luminosity of the SFRs in 
our sample (solid line) and that of a larger sample of SFRs 
(Palagi et al. 1993) (dashed line). Both histograms are nor- 
malized to the total number of sources. 



Fig. 3. The distribution of the integrated fluxes of the H2O 
masers in our sample as measured in the Feb. 2007 observ- 
ing session (solid line) , and that of a larger sample of SFRs 
(a merger of Comoretto et al. 119901 Brand et al. 119941 and 
Valdettaro et al. l200ip (dashed line). 



The available bandwidth varied between 3.125 MHz and 
25 MHz (at 22 GHz this corresponds to a velocity resolution 
between 0.041 km and 0.658 km s^^). In February 1997, 
the VLA-1 chip correlator was replaced by a new one based 
on the NFRA correlator chips and boards (Bos ll99T|) . The 
new autocorrelator has 2048 delay channels and 160 MHz 
maximum bandwidth in its standard configuration. In or- 
der to maintain compatibility with previous observations, 
we have continued to use only 1024 channels and a 10 MHz 
bandwidth, giving a velocity resolution of 0.132 km s~^. 
Finally, owing to extensive maintenance works on the ra- 
diotelescope, there is a gap in the observations from April 
1996 to February 1997. 

The telescope pointing model is typically updated a few 
times per year, and is quickly checked every few weeks by 
observing strong maser sources (e.g., W3 OH, Orion-KL, 
W49 N, Sgr B2 and W51). The pointing accuracy was al- 
ways better than 25"; from about mid-2004 the rms resid- 
uals from the pointing model were about 8"-10". 

Observations were taken in position-switching mode, 
with both ON and OFF scans of 5 min duration. The OFF 
position was taken 1.25° E of the source position to rescan 
the same path as the ON scan. Typically, one ON/OFF 
pair was taken, though for weak sources this was repeated 
several times. The typical 1 a noise-level in the spectra is 
-1.5 Jy. 



3.1. Calibration 

The calibration of the intensities and the velocities is one 
of the main concerns for variability studies. 

The antenna gain as a function of elevation is deter- 
mined by daily observations of the continuum source DR 21 
(for which we assume a flux density of 16.4 Jy after scaling 
the value of 17.04 Jy given by Ott et al. 1994 for the ratio of 
the source size to the Medicina beanfl) at a range of eleva- 
tions. Antenna temperatures were derived from total power 
measurements in position switching mode. Integration time 
in each position is 10 sec with 400 MHz bandwidth. The 
zenith system temperature is about 120 K in clear weather 
conditions. 

The daily gain curve was determined by fitting a poly- 
nomial curve to the data, which was then used to convert 
antenna temperature to flux density for all spectra taken 
that day. From the dispersion of the single measurements 
around the curve we find the typical calibration uncertainty 
to be 19%. On the few days for which no separate gain curve 
was measured, we applied that which was closest in time, 
and estimated a corresponding calibration uncertainty of 

^ In previous works and up to April 2003 the flux density of 
18.8 Jy was adopted for DR 21, following Dent ([T972|l . In this 
paper, all data were recalibrated using the value given by Ott 
et al. (p94l) . 
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7%. The average overall calibration uncertainty is therefore 
estimated to be ~20%. 

To reduce the effects of atmospheric attenuation vari- 
ations, observations were obtained in good weather con- 
ditions and at elevations > 30°, or around source merid- 
ian transit for low declination sources. For several sources 
we compared our intensities with those available in the 
literature and close in time with our observations, virtu- 
ally always finding very good agreement within the above- 
quoted uncertainty (e.g., Valdettaro et al. I2002p . A good 
agreement was also found comparing observations taken at 
Effelsberg in a program to study H2O masers associated 
with late-type stars which was run jointly by the Medicina 
and Effelsberg radiotelescopes. 

As for the calibration of the velocity scale and its sta- 
bility over time, checks were made with quasi simulta- 
neous Effelsberg observations showing a good agreement. 
However, it is difficult to quote a definite uncertainty con- 
sidering the variable nature of the sources and the not si- 
multaneous observations. To make sure that the velocity 
changes that will be reported are real effects and not due 
to instrumental changes (e.g., local oscillator drifts in our 
system) and to evaluate the uncertainty, we performed in- 
ternal consistency checks within the sources of our sample. 
The source G32.74-0.08 is characterized by a single, nar- 
row and intense velocity component. The observed velocity 
of the peak of the emission displays a maximum deviation 
from a linear fit over the entire period of < 0.1 km s~^, a 
value that will be used as maximum possible uncertainty 
on the velocity throughout our patrol. 



4. The database 

For the purposes of the present study, the standard reduc- 
tion procedures were applied to the observed spectra before 
including them in the database: 

(a) the edges of the available band were deleted to limit 
the spectra to the flatter part of the band. This means 
that the velocity coverage shown may be smaller (gen- 
erally by 20-30%) than the nominal bandwidth of the 
autocorrelator; 

(b) the spectra taken within 10 days were averaged, unless 
a significant difference was found; 

(c) a polynomial fit to the channels free of line emission was 
removed from the spectra; 

(d) rebinning to 0.3 km was applied. 

The whole set of Medicina 32-m observations is archived 
m CLAS^I format. The raw data can be accessed through 
contact with the authors. 

The amount of information present in the database is 
so large that a concise way of presentation is essential for 
a rapid visual appreciation of the properties of maser vari- 
ability. Below we describe our selection. All plots and fig- 
ures presented can be found in electronic form in the online 
appendix. 



^ CLASS is part of the GILDAS software de- 
veloped at IR AM and Observatoire de Grenoble 
| http://www.iram.fr/IRAMFR/GILDAS| ) 



4.1. The spectra 

The most obvious presentation is a time sequence of spec- 
tra. These are arranged in a compact form as page-plots 
in which each page contains 18 spectra, with time running 
from top to bottom and from left to right. In each spec- 
trum, the date of observation is reported in the top left 
corner, and the number of days elapsed since the first ob- 
servation of the source is shown in the top right corner. 
The velocity scale is the same for all spectra of each source 
and is bracketed by the minimum and maximum velocities 
at which emission has been detected (above the 5 a level) 
throughout the monitoring period. Given the large range of 
intensities present during our patrol, we present the spec- 
tra using an autoscaled flux density scale. These plots are 
convenient to view the full details of each spectrum except 
for those with bursts, although they have the disadvantage 
that the flux density scale changes from spectrum to spec- 
trum and the visual impression of the variability is lost. 

As an example we show in Fig. 2] the flrst page of the 
spectra for the source W75-N. 

The following two additional presentations of the spec- 
tra can be downloaded from our WEB pagetQ. 

(1) Fixed-linear flux density scale, where the maximum flux 
density is determined by the maximum value observed 
during our monitoring. These plots are convenient for 
a quick estimate of the variability, but compress the 
spectra too much when the source presents bursts with 
intensities orders of magnitude larger than the "quies- 
cent" value. 

(2) Same as in (1) but with logarithmic flux density scale. 
The minimum flux density is equal to the 5 a level in 
each spectrum. 

4.2. The velocity-time-flux density plots 

A more concise visualization of the time variation of the 
maser emission can be obtained by showing the flux density 
versus velocity and time in the same diagram. 

Since the observed dates are not distributed evenly in 
time, some kind of interpolation between two adjacent ob- 
servations must be used. As no assumption can be made 
on the evolution of the H2O maser components, we chose a 
linear interpolation. This procedure produces an apparent 
increase in the lifetime of a feature when long time-intervals 
between two consecutive observations occur. 

Two versions of these plots are presented: 

(1) full plots, with velocity coverage extending to the max- 
imum velocity covered during the observations; 

(2) zoomed plots, with velocity coverage limited to the 
part of the spectrum where emission has been detected 
(above the 5 a level) at least once during the patrol. 

In these plots, only emission above the 5 a level is 
shown. In both plots the vertical solid line indicates the 
velocity of the thermal molecular gas (either CO, CS or 
NH3). In the zoomed plots the additional dashed line repre- 
sents the mean velocity derived from the histogram of the 
rate-of-occurrence (see Sect. 4.5). The time-scale on the 
left is in days starting from the flrst observation. To con- 
vert this scale to real dates one should use the date given 

^ http://www.arcetri.ci8tro.it/~starform/water_maser_v2.html 
http://www.ira.inaf.it/papers/masers/water_maser_v2.html 



M. Felli et al.: Variability of H2O masers in SFR 



7 



BOS 

3 



> 1 1 1 1 1 1 


■ 1 ' 1 J 


■- 

1 1 1 1 1 1 1 1 1 


- 1 1 1 1 i 1 1 1 



D „ i'l-S£P-lMS 
BOD 



K)0 - 

quy - 

SOD - 



W75-N 



1"^ 



i36 




Ilia 



i«o - 



J 




ET — i 1 — I — 1 — r 



1 smpi 



2M 31-JIIW-l 
1M6 






1 1 1 


P w 
2i-IHH*-ia8a 


r' ■ 1 ' ' 

; j 


j 


2(5 
a-M(5V-''-flM J3J 






□ ;(> 

(H-FEH-:9eC V161 



1«0 



J 



i9-jt>j-i9ai 



ItCC! !^ 




1 — ' — I — ' — I — r 



T — I — 1 — I — r 



1 I 1 1 1 l'\ 1 . II 


c 

— = — 1 — 1 — 1 — 1— 


20 
T 1 — r 


1^ 




r7-JW-1Mi ITfiS 

C — ' — I — ) — I — ITT — I — I — I — I 



V,^R (km/s) 



Fig. 4. The first 18 spectra of the source W75-N with an autoscaled flux density scale. Time runs from top to bottom 
and from left to right. In each spectrum the actual date of observation is given above the top left corner, and the number 
of days elapsed since the first observation of the source above the top right corner. The velocity scale is the same for all 
spectra. 



in the first spectrum of the figure with the time sequence of As an example of the velocity-time-flux density full and 
spectra, e.g., Fig.|3]for W75-N or the date listed in Table 1. zoomed plots, in Fig.[5]-{B]we show those of the source W75- 
These plots are particularly useful to find possible velocity N. 
drifts and to separate steady components from bursts of 
short duration. 
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W75-N 




50 



Fig. 5. The velocity-time-flux density full plot of the source W75-N. The vertical solid line indicates the velocity of the 
associated thermal molecular gas. The flux density scale is shown by the bar on the right. In this bar the three lines 
give the flux density of the drawn contours. The area in black indicates that the corresponding velocity range was not 
observed. (This figure is in colour both in the electronic edition and in the online appendix.) 



4.3. The light curve of the maser emission 

A plot of the peak fiux density as a function of time is 
often used to study the temporal behaviour of the maser 
emission. However, the peak flux density is meaningful only 
if one narrow component is present throughout the entire 
observing period. It has less physical meaning for complex 
spectra with multiple components of comparable flux densi- 
ties varying in an uncorrelated fashion. In this case, the in- 
tegrated flux density can only be used with the understand- 
ing that this represents the contribution of all components 
and may be representative of the energy input from the 
YSOs. However, since each velocity component may vary 
independently, the relative contribution of each to the light 
curve is lost in this plot. 

An example of the light curve of the integrated flux 
density for the source W75-N is shown in Fig. [T] 

4.4. The upper and lower envelopes and the mean spectrum 

Another meaningful description of the degree of maser vari- 
ability is given by the comparison between the upper and 
lower envelopes and the mean spectra of the sources over 
the whole period of observation. 



W75-N 

4000 F~i ' ' ' I ' ' ' r- 




L-j I 1 I I I I 1 \ I 1 1 \ 

2000 4000 6000 

time (days) 



Fig. 7. The integral of the flux density over the observed 
velocity range as a function of time for the source W75-N. 



The upper envelope corresponds to the maximum flux 
density ever reached at each velocity and is obtained by 
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Fig. 6. The velocity-time-flux density zoomed plot of the source W75-N. The vertical solid line indicates the velocity of 
the associated thermal molecular gas. The vertical dashed line marks the mean velocity derived from the histogram of 
the rate-of-occurrence. The flux density scale is shown by the bar on the right. In this bar the three lines give the flux 
density of the drawn contours. The area in black indicates that the corresponding velocity range was not observed. (This 
figure is in colour both in the electronic edition and in the online appendix.) 



finding the maximum flux density in each velocity bin. 
Its integral and the derived luminosity, Lh2o(up), repre- 
sent the maximum emission that could be produced by the 
source if all the velocity components were to emit at their 
maximum level and at the same time. 

The mean spectrum represents the mean flux density at 
each velocity, and is obtained by computing the arithmetic 
mean of the flux densities in each velocity bin, assigning the 
same weight to the individual spectra. In this average the 
flux density is arbitrarily set to zero if it is below the 5 a 
noise level in the spectrum in question. Clearly, the steadier 
components dominate this spectrum. 

Finally, the lower envelope identifies the components 
that are always active and is obtained by finding the min- 
imum flux density in each velocity bin, again setting it to 
zero if it is below the 5 a noise level in each spectrum. 
When the lower envelope is zero over the whole velocity in- 
terval, it indicates that no velocity component is constantly 
present above a 5 cr level during the monitoring period. This 
does not necessarily imply that the maser is quiescent at 
all velocities at some epoch, since the emission may occur 
at different velocities at different times. 



An example of the upper and lower envelopes, and mean 
spectrum of the source W75-N is shown in Fig. [H 

4.5. The rate-of-occurrence plots 

To identify the velocity components which are more fre- 
quently present in the spectra (independently of their flux 
density), we plot the rate-of-occurrence of maser emission 
above the 5 a noise level as a function of velocity. To pro- 
duce this plot we have set a counter for each channel that 
increases by one unit every time the flux density in the 
channel is greater than the Scr noise level of the spectrum. 

An example of the rate-of-occurrence histogram of the 
source W75-N is shown in Fig. (5] 

Two histograms can be produced: 

(1) the percentage of detections with respect to the total 
number of observations (dotted line histogram, scale to 
the right); 

(2) the number of detections (solid hne histogram, scale to 
the left). 

The two histograms may differ slightly in velocity ranges 
where the time coverage is not equal for all channels, pro- 
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Fig. 8. The upper (red) and lower (blue) envelopes and 
the mean spectrum (green) of the source W75-N measured 
during our monitoring. The vertical solid line marks the ve- 
locity of the associated thermal molecular gas. The vertical 
dashed line marks the mean velocity derived from the his- 
togram of the rate-of-occurrence. (This figure is in colour 
both in the electronic edition and in the online appendix.) 
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Fig. 9. The rate-of-occurrence plot for the source W75- 
N. The scale to the right refers to the dotted histogram, 
the scale to the left to the solid line histogram. The verti- 
cal solid line marks the velocity of the associated thermal 
molecular gas. 



duced by using different band widths during our monitoring. 
From this histogram we can derive the mean velocity which 
is shown in Figs. [6] and [8l 



5. Conclusions 

We present the observational results of a systematic study 
extending for almost 20 years of the H2O maser variability 
in 43 SFRs with luminosities of the associated IR sources 
between 20 and 1.5x10^ L©. This database provides the 
backbone for the discussion of the main long-term prop- 
erties of maser emission that will be presented in a forth- 
coming paper. We have identified several ways to describe 
graphically the main aspects of the H2O maser emission. 
These include: 



(a) 



(b) 



(c) 



(d) 
(e) 



(f) 



the spectra in a compressed form, with an autoscaled 
flux density scales (the same sets of plots, but with fixed 
linear and logarithmic flux density scales, are also avail- 
able from our WEB pages); 

the velocity-time-flux density plots, which conveniently 
describe the morphology of the variability of the maser 
emission. These diagrams are particularly useful for rec- 
ognizing the presence of possible velocity drifts and sep- 
arating steady components from bursts of short dura- 
tion; 

the velocity-integrated flux density as a function of 
time. This light curve describes the overall emission of 
the maser components associated with a SFR; 
the mean spectrum; 

the upper and lower envelopes of the maser emission. 
The upper envelope represents the maximum emission 
that the SFR could produce if all the velocity compo- 
nents were present simultaneously and emitting at their 
maximum rate. Similarly, the lower envelope pinpoints 
the steady components and their lowest level (but > Scr) 
of emission; 

the number of times that the maser emission is well 
above the noise (> 5a) as a function of velocity 
is directly gauged by the histogram of the rate-of- 
occurrence. 
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